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ABSTRACT 



Results on steady TeV 7-ray point source search using data taken from the 
Tibet HD (Feb. 1997 - Sep. 1999) and Tibet III (Nov. 1999 - Oct. 2001) 
arrays are presented. From 0° to 60° in declination, significant excesses from the 
well-known steady source Crab Nebula and the high state of the flare type source 
Markarian 421 are observed. Because the levels of significance from other posi- 
tions are not sufficiently high, 90% confidence level upper limits on the flux are set 
assuming different power law spectra. To allow crosschecking, two independently 
developed analyses are used in this work. 

Subject headings: gamma-rays, all-sky survey, point sources 
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1. Introduction 

In the last 15 years, Imaging Atmospheric Cherenkov Telescopes (IACT) have discov- 
ered more than 10 sources in TeV energy range populating from pulsar nebula, supernova 
remnants, starburst galaxies, AGNs of BL Lac type and even a newly found unidentified 
source (Ong 2003; Horan et al. 2004; Weeks 2003; Volk 2003), and it is remarkable that 
AGNs make up a large percentage of these sources. To understand the mechanism of the 
TeV 7-ray emissions from diversity of objects, and more fundamentally, to understand the 
origin and acceleration of cosmic ray (CR), it is important that the number of TeV sources 
will be largely increased. However, to improve the statistics of TeV sources is not a trivial 
task for IACT due to its small field of view and low duty cycle, which particularly limit its 
power in detecting flare type sources and unknown sources. Therefore, a complementary de- 
tecting technology is necessary. Although their sensitivity for detecting dim 7-ray sources is 
lower than that of the IACT, Extensive Air Shower (EAS) experiments, such as the Tibet air 
shower array (Amenomori et al. 1999; Amenomori et al. 2000; Amenomori et al. 2003) and 
Milagrito/Milagro (Atkins et al. 1999; Atkins et al. 2003; Atkins et al. 2004) experiments, 
have resulted in successful observation of 7-ray emissions from standard candle Crab nebula 
and from transient sources such as Mrk 501 and Mrk 421. Their characteristic abilities in 
high duty cycle and large field of view allow them to simultaneously monitor a larger area in 
space over continuous time. Most importantly, the EAS technique can potentially provide 
useful information for IACT observatories for their further dedicated observations. Such 
follow-up observations have been performed by Whipple experiment (Falcone et al. 2003a; 
Walker et al. 2003b) guided by the hot spots seen from TeV 7-ray all-sky survey results pub- 
lished by the Tibet air shower array (Amenomori et al. 2001a; Cui et al. 2003) and Milagro 
(Sinnis et al. 2002) experiments. 

The satellite experiment EGRET on board the Compton Gamma Ray Observatory 
(CGRO) pioneered the all-sky 7-ray survey in the energy range of 20 MeV to 30 GeV 
(Hartman et al. 1999). The successful detection of the diffuse 7-ray emissions from the 
galactic plane and a peak near 70 MeV seen in a rather hard 7-ray spectrum supported the 
supernova acceleration model, according to which TeV 7-ray emissions were also expected. 
In addition, EGRET discovered 271 7-ray point sources (Hunter et al. 1997), among which 
about 100 were identified as SNRs, AGNs, etc., while the rest has remained unidentified. The 
diverse types of 7-ray emissions, which included galactic sources and extragalactic sources, 
greatly improved our knowledge of both CR physics and astrophysics. 

The all-sky survey remains one of the major concerns for ground-based observatories 
(Weeks et al. 1989; Alexandreas et al. 1991; Mckay et al. 1993; Wang et al. 2001; Aharonian 
et al. 2002a; Aharonian et al. 2002b; Atkins et al. 2004; Antoni et al. 2004). For example, 
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AIROBICC made a northern sky survey (Aharonian et al. 2002a) on 7-ray emissions for 
energy above 15 TeV; because there was no compelling evidence for a signal source, an 
absolute flux upper limit between 4.2 and 8.8 Crab was obtained for a declination (Decl.) of 
0° and 60°. With the IACT technique, within the TeV energy range, the HEGRA telescope 
was used to performe a survey (Aharonian et al. 2002b) in one quarter of the galactic plane 
in the longitude range from —2° to 85°. As a result of negative evidence, upper flux limits 
on each individual known point source located in this region were obtained between 0.15 
Crab and several Crab. Most recently, Milagro updated its results (Atkins et al. 2004) in 
the northern sky survey in TeV energy range, and pushed the average flux upper limit down 
to a level between 275 and 600 mcrab (4.8 x 10~ 12 to 10.5 x lO^cm-V 1 ) above 1 TeV 
for source Decl. between 5° and 70°. As an independent experiment, the Tibet air shower 
array has similar sensitivity, and covers almost the same energy ranges and field of view as 
those experiments do. The results from the Tibet air shower array should therefore provide 
important information for crosschecking and confirmation. 



2. Tibet Air Shower Observatory 

The Tibet air shower array has been conducted at Yangbajing (90.522°E, 30.102°N; 
4,300 m a.s.l) in Tibet, China since 1990. The Tibet I array (Amenomori et al. 1992), which 
consisted of 49 scintillation counters forming a 7 x 7 matrix of 15 m span, was gradually 
expanded to the Tibet II array occupying an area of 36,900 m 2 by increasing the number 
of counters from 1994 to 1996. Both Tibet I and II have the same mode energy (the most 
probable gamma ray energy observed from the standard candle Crab) of 10 TeV. In order to 
observe TeV CRs, in 1996, part of the Tibet II array covering an area of 5,175 m 2 was further 
upgraded to a high density (HD) array (Amenomori et al. 2001b) with 7.5 m span. Soon 
after decreasing the mode energy to 3 TeV, the HD array observed multi-TeV 7-ray signals 
from the Crab Nebula (Amenomori et al. 1999) and Mrk 501 (Amenomori et al. 2000). To 
increase the event rate, in 1999, the HD array was enlarged to cover the center part of the 
Tibet II array as the Tibet III array (Amenomori et al. 2001c). The area of Tibet III has 
reached to 22,050 m 2 (Figure 1). 

As can be shown in Figure 1, the Tibet HD and Tibet III have identical structure except 
the array size and shape. Each counter has a plastic scintillation plate (BICRON BC-408A) 
of 0.5 m 2 in area and 3 cm in thickness, and is equipped with a fast-timing (FT) photomul- 
tiplier tube (PMT; Hamamatsu H1161). A 0.5 cm thick lead plate is put on the top of each 
counter in order to increase the array sensitivity by converting 7-rays into electron-positron 
pairs in the shower (Bloomer et al. 1988; Amenomori et al. 1990). The angular resolution 
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of the HD and the Tibet III are about 0.9° in the energy region above 3 TeV, as estimated 
from full Monte Carlo (MC) simulation (Amenomori et al. 1990; Kasahara 2003) and ver- 
ified by the moon shadow measurement from observational data (Amenomori et al. 1993; 
Amenomori et al. 2001c; Amenomori et al. 2003). The trigger rate is about 105 Hz for the 
HD and 680 Hz for the Tibet III. In this work, the sample includes data obtained by running 
the HD array for 555.9 live days from Feb. 1997 to Sept. 1999 and the Tibet III array for 
456.8 days from Nov. 1999 to Oct. 2001. 

The event selection was done by imposing the following four criteria on the reconstructed 
data: i) each shower event should fire four or more Fast-Time (FT) detectors recording 1.25 
or more particles; ii) the estimated shower center location should be inside the array; iii) 
the sum of the number of particles per m 2 detected in each detector Pft should be larger 
than 15; iv) the zenith angle of the incident direction should be less than 40°. After applying 
these cuts and a data quality cut, about 40% of the shower events were selected, results in 
a total number of events to be about 7.0 x 10 9 . 



3. Analysis 

Because the Tibet air shower array can not distinguish a 7-ray induced shower event 
from the overwhelming CR background shower events, when tracing and counting the number 
of events in an "on-source window" centered at a candidate point source direction with a size 
at the level of angular resolution, the number of background events must be estimated from 
the observational data recorded in the side band, which is usually referred to as "off-source 
window" . 

Sitting on an almost horizontal plane, the Tibet HD and III have almost azimuth- 
independent efficiency in receiving the shower events for any given zenith angle. The equi- 
zenith angle method was therefore developed. In brief, simultaneously collected shower 
events in the same zenith angle belt can be used to construct the "off-source windows" and to 
estimate the background for a candidate point source located in the same zenith angle. This 
method can eliminate various detecting effects caused by instrumental and environmental 
variations, such as changes in pressure and temperature which are hard to be controlled and 
intend to introduce systematic error in measurement. In reality, the actual azimuth angle 
distribution deviates from a uniform one mainly due to the fact that the Tibet HD and III 
arrays are situated on a slope of about 1.3° from the southeast to the northwest. Furthermore, 
a geomagnetic effect in the northern hemisphere causes unequal efficiency between north 
and south (Ivanov et al. 1999). Together with other possible unknown effects, all the steady 
effects lead to an approximately 2.5% higher event rate from the southerly direction than 
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from the northerly direction. This non- uniformity must be properly accounted for in the 
following analyses. 

For the purpose of crosschecking, two independently developed analyses based on the 
equi-zenith angle method are used in this work. One is dedicated to a point source search, 
while the other in addition is sensitive to large-scale anisotropy of the intensity of the CR. A 
MC study was used to carefully compare the two methods, and the search sensitivities were 
found to be very comparable. However, the number of the background events estimated 
by the two methods is different due to the different ways used in choosing the off-source 
windows (see subsections below in detail). As a result, significance values calculated from 
the two methods differ by about 0.5 a statistically. 



3.1. Method I (short distance equi-zenith angle method) 

As a good approximation, the large-scale anisotropy of the CR intensity can be neglected 
if the "off-source windows" are chosen to be close enough in distance to the "on-source 
window", since the CR intensity only changes slowly in the sky. In this analysis, to make 
sure not to miss any possible unknown source, the surveyed sky has been over-sampled 
by the following way: the sky is divided into 0.1° x 0.1° cells, from 0° to 360° in Right 
Ascension (R.A.) and from 0° to 60° in Decl. Pointing to the center of those cells, cones 
with a half opening angle of 0.9° (for E > 3 TeV) or 0.4° (for E > 10 TeV) are tested as 
on-source windows. These cone sizes are chosen to maximize the source search sensitivity 
(Chen et al. 2004). 10 off-source windows of the same shape are symmetrically aligned 
on both sides of the on-source window, at a zenith angle-dependent step size in order to 
maintain an angular distance of 2° between the neighboring off-source windows. In order to 
avoid signal events related to an on-source window being wrongly counted as a background 
event, the two off-source windows closest in distance to the on-source window are required 
to be placed at twice of the step-size from the center of the on-source window. Denoting 
the number of events in on-source window as ./Von, and the number of events in the ith 
off-source window as AoFFi, without considering the effect from non- uniform azimuth angle 
distribution, the number of background events can be estimated as 

N BG = E '=i NoFFi (1) 
10 K ' 

When the zenith angle is less than 4.2°, the 10 off-source windows defined above become 
to be overlapped. Therefore, data obtained with a zenith angle less than 4.2° are not used. 
Since the azimuth angle distribution is not completely flat as discussed earlier, the averaged 
number of events from off-source windows, N BG , as calculated in equation (1) is not an exact 



description of the background number in the on-source window. Therefore Nbg needs to 
be corrected. Giving the fact that azimuth angle distribution remains stable over a period 
of time that is much longer than the time scale of a day, the effect of the tiny violation of 
the exact equi-zenith angle condition can thus be determined from the observational data 
which were taken only a few hours apart from the exposure of the on (off)-source window(s). 
As a matter of fact, the previous measurement is repeated 35 times (the relevant "on" and 
"off" source windows are hereafter called "dummy on" and "dummy off" source windows). 
Together with the on-source window, 35 dummy on-source windows on the orbit which 
has the same Decl. as the on-source window are defined at regular intervals along R.A. 
direction, i.e., the dummy on-source windows are located at ojon + 10° x j, where «on is 
the R.A. location of the on-source window, and j an index of the repeated measurements 
which runs from 1 to 35. N' ON j and N# G j represent the number of events in the dummy 
on-source window and its estimated background from the jth measurement, respectively. 
The correction factor, i], due to non-uniform azimuth angle distribution is given by 



As for the Tibet observation, 77 deviates from 1 at an order of 10 3 depending on the Decl. 
of the on-source window. 

Finally, the corrected estimation of the number of background events is obtained as 



It is clear that such a correction is free from the systematic uncertainty due to the slow 
variation of instrumental and environmental effects. Denoting the number of excess events 
in an on-source window as N$, its uncertainty, AiVs, is given by AiVs 2 = AiVoN 2 + AiVg G 2 , 
where AiV N 2 is well approximated by iVg G , and AiVg G is propagated from the statistic 
error of 77 and iV BG , according to Equation (3), AiV^ G 2 = N^ G 2 (^- + Considering 
the number of off-source window and dummy on(off)-source window, it is not difficult to 
calculate AiV* G as N^ G 2 [(± + ^X^) + (£) (j^)} = 0.1314iY^ G . Then we obtain AN S = 
1.0637^/iVg G . Because both iV Q N and N^ G are statistically large numbers, the significance 
value can be simply calculated as 



7] = 



(2) 



iV*BG = V^BG 



(3) 
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3.2. Method II (all distance equi-zenith angle method) 

Unlike the method described in the previous sub section, this method attempts to exploit 
the statistics as much as possible from two approaches. 

First, shapes of non-uniformly distributed azimuth angle distribution from whole data 
set are used to do the azimuth correction. To be more precise, after event selection, for 
each 1° zenith angle interval between 0° to 40°, azimuth angle distributions are filled and 
normalized, so that the averaged function value is one. Inversing the function and using it 
as an event weight would keep the total number of events unchanged but make the azimuth 
angle distribution flat. 

Second, all events in the equi-zenith angle belt except those inside the on-source window 
are taken as off-source events. In this case the large-scale anisotropy effect can no longer be 
ignored, which leads to the attempt to fit simultaneously the relative CR intensity over the 
sky within the detector's field of view. 

The idea of this method is that at any moment, for all directions, if we scale down (or 
up) the number of observed events by dividing them with their relative CR intensity, then 
statistically, those scaled observed number of events in a zenith angle belt should be equal 
anywhere. A x 2 function can be built accordingly and the relative intensity of CRs in each 
direction can be solved by minimizing the x 2 function. 

As in the previous method, the celestial space from 0° to 360° in R.A. and from 0° to 60° 
in Decl. are binned into cells with a bin size of 0.1° in both the R.A. and Decl. directions. 
In the observer's coordinates, the zenith angle 9 is divided from 0° to 40° by a step size 
of 0.08°, and the azimuth angle is binned by a zenith angle dependent bin width (0.08°/ 
sin(0)). For every local sidereal time (LST) interval bin (24s) m, a cell in (9, 0) space (n, /) 
is mapped to a celestial cell (««, 5j), i.e. R.A. bin i and Decl. bin j, through two discrete 
coordinate transformation functions t(m,n,l) and j(m,n,l). Therefore, at a certain LST 
bin m, the number of events accumulated in zenith angle bin n and azimuth angle bin I is 
directly related to the CR intensity in cell of celestial space. 

Denoting the number of observed events after azimuth angle correction as ./Vobs(^, n, I), 
and the relative intensity of CR as /(i,j), the equi-zenith angle condition leads to the following 
X 2 function: 

2 _ [iVoBsKn, /)//(*, j) - Ei>*i(NoBs(m,n,l')/I(i', f)) / Ev^lf 

X m,n,l AW™, / ^ j) + £, ¥ , (^OBsfo 71, V) / P(l>, j')) / (£*'_,/ l) 2 

(5) 
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Here, is mapped from (m,n,l) by the above-mentioned transformation functions. 

From equation (5), j) and its error AI(i, j) can be solved numerically by the iteration 
method. As the dimensionless variable is only a relative quantity and because the 

detecting efficiency along the Decl. direction can not be absolutely calibrated, we apply the 
constraint that averaged in each Decl. slices is equal to one. 

On the other hand, N(i,j), the number of observed events in the celestial cell is sim- 
ply a summation of the relevant number of events counted in local coordinates A^ Bs( m ; n i 0> 
i.e.,N(i,j)= Yl N OBS (m,n,l). 



As the relative intensity of the CR contains the contribution from large-scale anisotropy 
of the CR, to remove this effect and to isolate the contribution from point sources, the 
measured CR intensity in a 10° belt along Decl. is projected to the R.A. direction, in 8° bin 
size. The center of the belts moves from 0° to 60°, in 1° step size. By smoothing those curves 
and subtracting them from the corresponding belt of 1° width, we obtain a CR intensity, 
icoRR(*j.7')> that is corrected for the anisotropy effect. 

With this corrected CR intensity, the number of excess events and their uncertainties 
in cell can be calculated as 



Taking into account the array's angular resolution, events are summed up from a cone 
with an axis pointing to the source direction, and the half opening angle is set as 0.9° (for 
E > 3 TeV) or 0.4° (for E > 10 TeV). All celestial cells with their centers located inside the 
cone contribute to the number of events as well as its uncertainty. Finally, the significance 
for an on-source window centered at cell can be calculated as 



?(m,n,/)=i 
3(m,n,l)=j 





AN s (i,j) = A/ C orr(«, j) • N(i,j)/I C oRR(i,j) 



(7) 



S(i,j) 




(IcORR(i',f) ~ 1) • N(l',f)/IcORR(t',f) 



(8) 
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4. Results 

The significance distributions from all directions are shown in Figure 2a and Figure 
2b for the two analysis methods, respectively. In the case of method II, the large-scale 
anisotropy must be measured and subtracted, the systematic uncertainty to the significance 
value due to this subtraction procedure is estimated to be 0.2 a by changing the bin size and 
the smoothing parameters. Figure 3 shows the intensity map containing the contribution 
from the large-scale anisotropy with a demonstration of its subtraction. 

The excellent agreement on the negative side with a normal distribution (Figure 2a, 
Figure 2b) indicates that the systematic effects are well under control for both analyses. As 
for the positive side, a wider shoulder exists with significance values greater than 4.0 a. The 
dominant contributions are due to the stable 7-ray source Crab and transient source Mrk 421 
which was in an active state from Feb. 2000 through Oct. 2001 (Amenomori et al. 2003). 
After removing their contributions, in such a manner that those cells that are less than 2° 
distance from the two sources are excluded, the dot-dashed histograms in Figure 2a and 
Figure 2b which show the significance distribution from the rest of the cells agree much 
better with a normal distribution. For reference, the Table 1 lists the directions of local 
maximal significance points with significance values greater than 4.5 a from either of the 
two methods, and those prominent directions contribute to the remaining small shoulder in 
Figure 2a and Figure 2b. 

Each significance value listed in Table 1 is calculated by equation (4) and (8) which is 
supposed to be used for a predefined on-source window. However, none of our prominent 
directions comes from a predefined one, as described in section 3.1, a large number of over- 
lapping on-source windows have been tried. In comparison with finding a high significance 
value from one predefined direction (e.g., Crab or Mrk421), such a procedure makes the 
probability being larger by a factor of about the number of trials, and therefore making the 
excess of on-source event number less significant. It should be mentioned that it is very 
difficult to count the number of trials exactly. To have order of magnitude estimation on 
the lower boundary of the number of trials, we can divide the solid angle of surveyed sky by 
the solid angle of one on source window, and get a number of about 7000 which still under 
estimates the number of trials. As an example, for the first prominent point in Table. 1, the 
significance value is 5.3 a from method 1(4.8 a from method II), enlarge the probability by 
7000 to account the number of trials effect end up with only 3.3(2.5) a in significance. To 
distinguish this significance from the previously calculated one, we call the former calculated 
one the pre-trials significance {S pre -t r iais) ■ This number of trials effect can also be understood 
from another point of view as the following: not including the Crab and Mrk421, both meth- 
ods observe four directions which have S pre ^ t riais greater than 4.5a. From pure background 
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MC simulation, the probability to observe four or more such directions is equivalent to a 
fluctuation of only 1.4(7. In conclusion, after taking into account the number of trials effect, 
the existence of the high significance point in Table 1 is well consistent with the background 
fluctuation. 

As can be seen in Table 1, two established sources Crab and Mrk 421 are detected 
in an angular distance not larger than 0.5° from their nominal positions, in agree with 
the 0.4° positional uncertainty estimated for the both methods. One of our previously 
reported (Cui et al. 2003) hot point with pre-trials significance value of 4.0 a located at 
(304.15°, 36.45°) is not listed in Table 1, because its pre-trials significance value is not im- 
proved after we further include more data samples into this analysis which increase the 
number of events by 60%: 4.1 a from Method I and 4.0 a from Method II. The interesting 
thing about this point is that it is close to the Cygnus arm and was found to be close to one 
of Milargo hot points (Walker et al. 2004), though the large angular distance of 1.5° indi- 
cates that the coincidence is not in favour of a point source hypothesis. Conclusive results 
will rely on further more observation. As a cross-check, Milagro hot point at (79.9°, 26.8°) 
(Atkins et al. 2004) is investigated in this work: the maximum pre-trials significant point is 
found at (79.3°, 25.9°) with 1.3 a from Method I, and (79.9°, 26.3°) with 1.7 a from Method 
II. In another word, no strong confirmation was obtained. On the other hand, except the 
agreement measurement with comparable sensitivity on Crab and Mrk 421 from both obser- 
vatories, no significant prominent direction from this work matches with Milagro hot points 
(Atkins et al. 2004) within the allowed source positional uncertainties. When not consider- 
ing the Crab and Mrk 421, we think the marginal coincidence and refutation of other hot 
points between Tibet Air Shower and MILAGRO could be attributed to the limited statistics 
and sensitivity of the two experiments. Whether these hot points are just due to statistical 
fluctuation or due to new TeV gamma sources will become clearer with the future improved 
statistic of observational data. 

As described previously, the significances from all other directions other than Crab and 
Mrk 421 are not high enough to definitely claim any existence of new TeV 7-ray point source, 
we set a 90% confidence level (CL) upper flux limit for all directions in the sky, except at 
the positions of Crab and Mrk 421. The upper limit on the number of events at the 90% CL 
is calculated firstly according to the number of excess events and its uncertainty for each cell 
following the statistic method given by Helene (Helene 1983). Then the effective detection 
area of Tibet Air Shower array is evaluated by full MC simulation assuming a Crab-like 
7-ray spectrum E~P with (3 = 2.6 for a set of Decl. values (0.0°, 10.0°, 20.0°, 30.0°, 40.0°, 
50.0°, 60.0°) and interpolated to other Decl. values between 0° and 60°. Taking into account 
the live time from observational data taking, we derive the 90% CL flux upper limit. Figure 
4 shows the maps of flux upper limits at 90% CL for energy E > 3 TeV and for E > 10 
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TeV (from Method I). Figure 5 shows the average flux upper limit along the right ascension 
direction as a function of Decl., which varies between 1.5 — 3.2 x 10^ 12 cm _2 s~ 1 for E > 3 
TeV and 0.2 — 0.4 x 10 _12 cm _2 s _1 for E > 10 TeV. Current integrated limits for energy 
above 10 TeV are the world-best ones which improve the results obtained from AIROBICC 
(Aharonian et al. 2002a) for energy above 15 TeV. For energy above 3 TeV, limits from this 
work are comparable with those obtained by Milagro (Atkins et al. 2004) for energy above 1 
TeV. The same procedure is applied to the cases of other power law indices for energy above 
3 TeV and above 10 TeV . The corresponding average flux limit for the other indices can 
be found in Figure 6 (from method II). Depending on the Decl. and power law index of the 
candidate source, the integral flux limits lie within 1.5 — 3.8 x 10 _12 cm~ 2 s~ 1 (E > 3 TeV) 
and 0.2 - 0.48 x lO^cm-V 1 (E > 10 TeV). 



5. Conclusion 

A northern sky survey for the TeV 7-ray point sources in a Decl. band between 0° to 60° 
was performed using data of the Tibet HD and the Tibet III air shower array obtained from 
1997 to 2001 with two independently developed analysis methods. The established Crab 
and Mrk 421 were observed. This indicates that the Tibet air shower array is a sensitive 
apparatus for TeV 7-ray astronomy and has a potential to find flare type 7-ray sources, such 
as BL-Lac type AGNs. In addition, the 6 other prominent directions were selected. However, 
more data will be needed before we can confirm or rule out these candidates. 

With the exception of Crab and Mrk 421, 90% CL flux upper limits are obtained from 
the rest of the positions under the hypothesis that a candidate point source are in power law 
spectra, with indices varying from 2.0 to 3.0. The integral flux limits lie within 1.5 — 3.8 x 
lO^cm-V 1 (E > 3 TeV) and 0.2 - 0.48 x lO^cm'V 1 (E > 10 TeV) depending on the 
Decl. and power law index of the candidate source. 

This work is supported in part by Grants-in-Aid for Scientific Research on Priority Areas 
(712) (MEXT) and by Scientific Research (JSPS) in Japan, and by a grant for International 
Science Research from the Committee of the Natural Science Foundation and the Chinese 
Academy of Sciences in China. 
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Table 1: Prominent directions with at least one pre-trials significance value greater than 
4.5 a from either of two methods. Each pre-trials significance values are calculated by the 
equation (4) or (8). 



No. 


Method 


R.A.(°) 


Decl.(°) 


N N 


NoFF 


N s 


AiV s 


Spre— trials 




I 


69.95* 


12.05 


828557.0 


823428.5 


5128.5 


965.2 


5.3 


1 


II 


70.25 


11.95 


825589.0 


821232.4 


4356.6 


906.2 


4.8 




I 


70.45 


18.15 


997441.0 


992485.7 


4955.3 


1059.7 


4.7 


2 


II 


70.45 


18.05 


994477.0 


990343.3 


4133.7 


995.2 


4.2 




I 


83.55 


22.25 


1087857.0 


1082101.0 


5756.0 


1106.5 


5.2 


3 


II 


83.35 


21.85 


1079728.0 


1074562.6 


5165.4 


1036.6 


5.0 




I 


88.85 


30.25 


1013683.0 


1008916.9 


4766.1 


1068.4 


4.5 


4 


II 


88.85 


30.25 


1185535.0 


1179647.4 


5887.6 


1086.1 


5.4 




I 


165.65 


38.45 


1170565.0 


1164237.2 


6327.8 


1147.7 


5.5 


5 


II 


165.55 


38.45 


1170632.0 


1164908.8 


5723.2 


1079.3 


5.3 




I 


221.65 


32.85 


1062942.0 


1057488.1 


5453.9 


1093.9 


5.0 


6 


II 


221.75 


32.75 


1191295.0 


1186765.1 


4529.9 


1089.4 


4.2 




I 


286.65 


5.45 


608892.0 


604934.2 


3957.8 


827.3 


4.8 


7 


II 


286.65 


5.55 


612041.0 


608322.9 


3718.1 


780.0 


4.8 




I 


309.85 


39.65 


1146305.0 


1141238.9 


5066.1 


1136.3 


4.5 


8 


II 


309.95 


39.65 


1146636.0 


1141544.7 


5091.3 


1068.4 


4.8 



*The last digital number in column of R.A. and Decl. are due to the way we divide the bin in the analyses. 

Note. — The meaning of the columns is (from left to right): sequence of prominent direction, analysis 
method, right ascension (J2000), declination (J2000), number of measured events in ON source window 
(Non), number of background events (Noff), event number excess in ON source window (Ns = N<jn — 
Noff), uncertainty on the event number excess (AN S ), and pre-trials significance S pre - t riais of deviation 
of Nqn from Noff- 
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Tibet III Air Shower Array (1999) 
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Tibet III In Area 
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■ FT Detector W/D-PM1 
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Fig. 1. — Schematical view of the Tibet III array operating at Yangbajing. Open squares: 
FT-detectors; filled squares: FT-detectors with a wide dynamic range PMT; filled circles: 
density detectors with a wide dynamic range PMT. We select air shower events whose cores 
are located within the detector matrix enclosed with the dashed line. 



Significance Distribution (Method I) 



Significance Distribution (Method II) 




Fig. 2. — The pre-trials significance distribution of all directions on the sky map from (a) 
Method I and (b) Method II. It should be mentioned that not all directions are statistically 
independent as the bin size is smaller than the angular resolution. The solid line is derived 
from all cells defined in analyses. The dot-dashed histogram with shading excludes those 
cells which have a distance to Crab or Mrk 421 shorter than 2°. The dashed histogram 
represents the best Gaussian fit to the data. 
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Map of Intensity | 




Fig. 3. — Upper panel: Map of the relative intensity of CR in the surveyed northern sky. 
Lower panel: A demonstration of the subtraction of the large-scale anisotropy of the CR. The 
anisotropy distribution in the narrow belt centered at 9° in Decl. is estimated by projecting 
the intensity measurement from a wider belt (a 10° width is chosen) in the R.A. direction 
with a bin size of 8°. This distribution is further smoothed and then subtracted from the 
narrow belt. 
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Fig. 4. — Maps of 90% CL flux upper limit assuming an energy spectrum of E 2 6 (from 
Method I). Upper panel: above 3 TeV; lower panel: above 10 TeV. 



- 20 - 



90% C.L Upper Limit (E>3TeV,E>10TeV) 
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Fig. 5. — R.A. direction averaged 90% CL upper limit on the integral flux above 3 TeV and 
10 TeV for a Crab like point source, i.e., with an energy spectrum of E~ 2 6 . The solid lines 
are obtained by Method I and the dashed lines by Method II. 
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90% C.L Upper Limit 
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Fig. 6. — R.A. direction averaged 90% CL upper limit on the integral flux above 3 TeV and 
10 TeV for different indices of power law spectra. As an example, results from Method II 
are shown in this plot. 



